Three-qubit Grover's algorithm using superconducting quantum interference devices 

in cavity-QED 
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We present a scheme for the implementation of three-qubit Grover's algorithm using four-level 
superconducting quantum interference devices (SQUIDs) coupled to a superconducting resonator. 
The scheme is based on resonant, off-resonant interaction of the cavity field with SQUIDs and the 
application of classical microwave pulses. We show that adjustment of SQUID level spacings during 
the gate operations, adiabatic passage, and second-order detuning are not required that leads to 
faster implementation. We also show that the marked state can be searched with high fidelity even 
in the presence of level decay and cavity dissipation. 

PACS numbers: 03.67.Lx, 42.50.Dv, 85.25.Cp 



I. INTRODUCTION 

The quantum algorithms work much more efficiently 
than their classical counter parts due to quantum super- 
position and quantum interference. For example, con- 
sider the search of an item in an unsorted database 
containing N elements. Classical computation requires 
0{N) steps to carry out the search. However, Grover 
showed that search can be carried out with only 0{^/N) 
steps P, 0. Thus, Grover's algorithm represents a 
quadratic advantage over its classical counterparts. 

Grover's algorithm has been realized using many phys- 
ical systems like NMR [3] , superconducting qubits [3, [1] 
and atom cavity QED systems Superconducting 
qubit cavity QED is an attractive approach for quan- 
tum information processing due to their strong coupling 
limit in microwave cavity as compared to atoms in cavity 
QED [Hi- SQUIDs have attracted much attention 
among the superconducting qubits, due to their long de- 
coherence time 12-14], design flexibility, large-scale in- 
tegration, and compatibility to conventional electronics 

The goal of this work is to implement three-qubit 
Grover's algorithm using four-level SQUIDs in cavity- 
QED. We consider a three-qubit phase gate, that reduces 
the number of quantum gates typically required for the 
realization of Grover's algorithm. Three-qubit Grover's 
algorithm is probabilistic [8J, as compared to two-qubit 
Grover's algorithm. Therefore, to achieve high success 
probability, we have to implement basic searching itera- 
tion several times. Implementation of three-qubit Grover 
search is much more complex as compared to two-qubit 
case. In our scheme, two lowest energy levels |0) and |1) 
of each SQUID represent logical states. The scheme is 
based on resonant, off-resonant interaction of cavity field 
with 1 2) — >■ 1 3) transition of SQUID and application of 
resonant microwave pulses. 

Our scheme does not require adjustment of SQUID 
level spacing during the implementation of Grover's 



search iteration, thus, decoherence caused by the adjust- 
ment of level spacing is suppressed. We do not require 
identical coupling constants of each SQUID with the res- 
onator and direct coupling between the level |1) and |0) 
[l7| . Grover's iteration time becomes faster due to reso- 
nant and off-resonant interactions as compared to second 
order detuning or adiabatic passage. 

Grover's iteration based on conventional gate decom- 
position consists of ten two-qubit phase gates [18| . How- 
ever, these ten two-qubit phase gates may be replaced by 
two three-qubit phase gates in our scheme which signif- 
icantly reduces the number of steps and the probability 
of error during the implementation. We also consider the 
effect of spontaneous decay rate from intermediate level 
1 3) and decay of cavity field during the implementation 
of Grover's iterations. 



II. GROVER'S ALGORITHM 

The basic idea of Grover's algorithm is as follows; we 
prepare input basis states in superposition state — 

(l/VlV) X^il^o^ N) by applying Walsh-Hadamard trans- 
formation. First we, invert phase of the desired basis 
state through unitary operator (called Oracle) and then 
invert all the basis states about the average. We con- 
sider the implementation of Grover's algorithm in terms 
of quantum logic networks as shown in Fig. 1. Any quan- 
tum logical networks can be constructed using quantum 
phase gates and single-qubit quantum gates. The single- 
qubit quantum gate for jth qubit can be written in Dirac 
notation as 



Hi^\ = cos ei 



^e{e-^nO){l\+ze^ni){0\). (1) 



For 6 — 7r/2 and <p = 0, we have H, 



Here 



ax is the Pauli rotation matrix whose function is to flip 
the state of qubit such that |0) -J> |1) and |1) |0) . 
For e = 7r/4 and ip = -7r/2, we have h[^^ = H^^^ 
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the sign of target state l^i, 92, 93)- The operator / — 
2 l^i, 92: 93) 92j 93I performs the unitary transforma- 
tion which can be implemented using three-qubit phase 
gate and single-qubit gate H2 — —i as shown 
in Fig. [T] The sign change operators for eight possible 
target states are given by 
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FIG. 1: Circuit diagram for three-qubit Grover's algorithm. 



which transforms each qubit into superposition state i.e., 
|0)^(|0) + |l))/^^and|l)^(|l)-|0))/^/2. 

The transformation for three-qubit quantum controlled 
phase gate can be expressed by 



— p(''V'55i,l'5g2. 1*93.1) 



91, 92, 93) 



(2) 



where jg)^, 1(7)2, and \q)^ stand for basis |0) or |1) of the 
qubit and (5gi,i, iJga.i, and ^^3,1 are the Kroneker delta 
functions. Thus, three-qubit quantum phase gate induces 
a phase e^^ only when all three input qubit are in state 
Three-qubit quantum phase gate operator for = tt 
can be written in Dirac notation as 



= |ooo) (ooo| -I- |ooi) (ooi| -f |oio) (oio| 

|011) (011| -I- 1100) (100| + |101) (101| 
-K|110)(110|-|lll)(lll|. 



(3) 



The three-qubit controlled phase gate can be used in- 
stead of involving series of two-qubit gates. This method 
not only simplifies the implementation but also reduces 
the probability of error. Figure. [1] shows the circuit di- 
agram of three-qubit Grover's algorithm based on three- 
qubit phase gate and two-qubit gates [19]. Consider that 
the initial state of three qubits is |000). Grover's algo- 
rithm can be carried out using the following three steps: 

Part 1 (W): Apply Walsh-Hadamard transformation 
W = Hi ® Hi ® Hf on each qubit. The resultant state 
is given by 



IV') = ^(|ooo) 



|100) 



|001) + |010) + |011) 
|101) + |110) + 



(4) 



Part 2 (C): In this step, consider the unitary op- 



erator C 



Co 



(called Oracle) which changes 



Now Oracle applies one of Cq^^q^^q^ operators on state 
given in Eq. Q and changes the sign of target state. For 
example, our target state is |001), then by applying Cqoi 
on state (|4]), we obtain the change of phase on target 
state |001) i.e.. 



IV'i) = ^(1000) - 1001) + 1010) + 1011) 
+ |100)-F|101) + |110) + |111)). 



(6) 



Part 3 (N): In this step, our goal is to find out the 
marked state |001). This can be accomplished through 
inversion about mean using the operator N = I — 
2\il])^{ip\. It is clear from Fig. [T] that the operator N 
can be written in terms of single-qubit gate and three 
qubit quantum phase gate 

N = {{-Hl®Hl®Hl))Q^{{Hl)~^®{Hl)-^®{Hl)-^}. 

(7) 

The combined operator G = —NCq^,q^,q^ is called 
Grover's operator. When G is applied to initial state 1-0) , 
k (« 7r\/]V/4) times, then the probability of searching 
target state becomes maximum [2|. 



III. BASIC THEORY 

Here, we consider rf-SQUIDs as a qubit that consists of 
a single Josephson junction enclosed by superconducting 
loop. The corresponding Hamiltonian is given by pOI] 



2L 



^ ,2Tr(j). , . 

Ejcos{-^), (8) 

00 



where G and L are junction capacitance and loop in- 
ductance, respectively. The conjugate variables of the 
system are magnetic flux </> threading the ring and total 
charge Qc on capacitor. The static external flux applied 
to the ring is (/)x and Ej = is the Josephson coupling 
energy. Here, Ic is critical current of Josephson junction 
and 00 = 5^ is the flux quantum. 
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FIG. 2: Resonant interaction of SQUID (1) and (2) with cav- 
ity field and off resonant interaction of target SQUID (3). 



The SQUIDs are biased properly to achieve desired 
four-level structure as shown in Fig. [5] by varying the 
external flux f2lj . The single- mode of the cavity field is 
resonant with |2) o |3) transition of SQUIDs 1 and 2. 
The evaluation of initial state |3) |0)^ and |2) u nder 
the effect of Hamiltonian dSI) can be written as [231 



|3) |0>, 
|2>|1>, 



cos(g,;i)|3)|0), 
cos{ga) \2)\1)^ 



isin{g,t) \2)\1)^, 
isin{g,t) \3}\0}^, (9) 



where |0)^ (|l)c) vacuum (single photon) states of 
the cavity field and gi {i = 1,2) is the coupling constant 
between the cavity field and |2) o |3) transition of the 
SQUID 1 and 2. 

The cavity field interacts off-resonantly with |2) ^ |3) 
transition of SQUID 3 (i.e., Ac = Wc — 1^32 >> 53) as 
shown in Fig. [21 Here, Ac is the detuning between 
1 2) o 1 3) transition frequency 1032 of SQUID 3, Wc is 
the frequency of resonator and g^ is coupling constant 
between resonator mode and |2) o |3) transition. In the 
presence of single photon inside the cavity, the evolution 
of initial state |2) and |3) are given by (HHI 



|2)|1),^ ^ e*^3*/^M2)|l)^, 



|3)|1), 



ii3l 



|3)|1) 



(10) 
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It is clear that phase shifts and e~*^ are induced to 
the states |2) and |3) of SQUID 3, respectively. However, 
states 1 2) |0)^ and |3) |0)^ remain unchanged. 

A classical microwave pulse resonant with \i) o |j) 
= 0,1,2,3) is applied to each SQUID. The evaluation 
of states can be written as [22] 



\i) — > cos{D,ijt) \i) — ie sin{flijt) \j) , 
|j) — )■ cos{nijt) \j) — ie^'^ sin{ilijt) \i) , 



where f2ij is the Rabi frequency between two levels |i) 
and \ j) and (f is the phase associated with classical field. 
It may be mentioned that resonant interaction between 
pulse and SQUID can be carried out in very short time 
by increasing the Rabi frequency of pulse i.e., amplitude 
of the pulse. 



IV. IMPLEMENTATION OF GROVER'S 
ALGORITHM 

Three SQUIDs shown in Fig. [2]are initially prepared in 
state |000). For notation convenience, we denote ground 
level as |1) and first excited state as |0) for SQUID 3 as 
shown in Fig. [2] 

Part 1 (W): To accomplish part loi Grover's algo- 
rithm, we apply single-qubit gate H{ to each SQUID as 
shown in Fig. [T] The single-qubit gate H{ is carried out 
through two-step process that involves an auxiliary level 
1 3) via method described in Ref [25|. We need two mi- 
crowave pulses of different frequencies. One is resonant 
to |1) O 1 3) transition frequency and other is resonant 
to |0) O 1 3) transition frequency. The desired arbitrary 
single-qubit gate (Eq. (1)) can be achieved by choosing 
the proper interaction time (i.e., 9 — r2o3i ) and phase 
If of classical microwave pulse resonance to |0) |3) 
transition. We achieve Hi by choosing fios^ = 7r/4 and 
ip — — 7r/4, as a result we obtained state given by Eq. 

©• 



Part 2 (C): In order to implement C = C, 



apply 

single-qubit gate H2 = —i cr^j by choosing ftost = 7r/2 
and 1^ = as shown in Fig. [1] Then apply three-qubit 
quantum controlled phase gate. The procedure for real- 
izing the three-qubit controlled phase gate is described 
as follows: 

Initially, cavity is in a vacuum state |0)^ and levels |2) 
and 1 3) of each SQUID are not occupied. 

Step 1. Apply microwave pulse (with ^isti = ^ and 
phase ip — TT, where ti is pulse duration) resonant to 
|1) |3) transition of the SQUID 1 to occupy level |3);^ . 
Cavity field interacts resonantly to the |2)-^ o |3)j^ tran- 
sition of SQUID 1 for time interval ti = TT/2gi such that 
the transformation |3);^ |0)^ — t- — i|2)^|l)^ is obtained. 
The overall step can be written as |1)^ |0)^ |2)^ |1)^. 
However, the state |0)j^ |0)^ remains unchanged. 

Step 2. Apply microwave pulse (with $^20^2 = t^/2 
and phase ip = 7r/2) to the SQUID 1 while a microwave 
pulse (with 1^20^2 — tt/2 and phase ip = —tt/2) to the 
SQUID 2. The transformation 12)^ {\0)^) \0)^ (- |2) J 
for SQUID 1 while |0)2 ([2)2) ^ [2)2 (- 10)2) for SQUID 
2 is obtained. 

Step 3. After the above operation, when cavity is in 
a single photon state |1)^, only the level |2) of SQUID 
2 is populated. The cavity field interacts resonantly to 
12)2 — )■ 13)2 transition of SQUID 2 for time t^ ~ 7r/2g2- 



We obtain transformation [2)2 |1) 



13)2 |0)^, while 



(11) 



12)2 |0) 

microwave pulse (with $^30^3 



„ ,1)210), and, 11)2 ID 



remain unchanged. Apply 
- tt/2 and phase p> = n) to 
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SQUID 2 to transform state |3)2 to i |0)2- The overall 
step can be written as 12)^ — 7> [0)2 |0)^. However, 
states 11)2 |0)^, 11)2 and [2)2 |0)^ remain unchanged. 
The evolution of the system after above three steps is 
given by 



|100)|0), 
|101)|0), 
|110)|0), 
1111) |0), 
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|200) 
|201) 
|210) 
1211) 



|020) 
1021) 
1010) 
1011) 



1000) |0), 

1001) |0), 

1010) |1), 

1011) |1)„ 



(12) 

It must be noted here that we have shown the evolu- 
tion of only four initial possible states out of eight states 
(See Eq. (4)) since the evolution of other four states i.e., 
|000) |0)^ , |001) |0)^ , |010) |0)^ , and |011) |0)^ is triviaL 

Step 4. Apply microwave pulse with 5112^4 = 7r/2 and 
phase if = ~tt/2 to SQUID (3) to obtain transformation 

11) 3 — ?► 12)3. After the above operation only level |2) of 
SQUID (3) is populated, when cavity is in single photon 
state. Now the cavity field interacts off-resonantly to 

12) 3 — > 13)3 transition of SQUID 3. It is clear from Eq. 
(7rAc)/.g|, the state [2)3 |1)^ changes 



10)3 |l)e and 



: n/2 
state 



pOj) that for 

to — 12)3 |1)^. However, the states 10)3 |0)^, 
12)3 |0)^ remain unchanged. 

Step 5. Apply microwave pulse (with fli2t4 - 
and phase (p = 7i'/2) to SQUID (3), as a result 
transformation [2)3 — > [1)3 is obtained. 

Step 6. Perform reverse of the operation mentioned in 
step 3. Apply microwave pulse (with ftsots = i'/2 and 
phase If — TT ) to SQUID 2. Wait for the time ^3 given 
in step 3, during which cavity field interacts resonantly 
to 12)2 — >■ 13)2 transition of the SQUID 2. Over all trans- 
formation can easily be written as |0)2 |0)^ — > |2)2 11)^. 
However, the states [1)2 |0)^, [1)2 |l)g and [2)2 |0)^ remain 
unchanged. The evolution of the system after applying 
steps 4-6 is given by 



|200) |1) 
1201) |1) 

1210) |1) 

1211) 11) 



|100) |0), 
1101) |0), 
1110) |0), 
-|111>|0), 



(14) 



After the application of three-qubit phase gate, the 
state given by Eq. Q evolves to 



IV'3 



1 

^ 2^72 
|100) 



(|000) + |001) + |010) + |011) 
101) + |110) - |111)). 



(15) 
This 



Then apply single-qubit gate i.e., H2 — —i Cx 
completes the part 2 (C) of the implementation. 

Part 3 (N): In order to implement operator N, apply 
single-qubit gate H^^, then apply three-qubit quantum 
controlled phase gate by repeating the above 8 steps. 
Finally, apply single-qubit gate Hi. It is clear that, 
Grover's operator G = —NCq^^q.^^q^ for eight objects 
{N = 8) can be implemented using four-level SQUIDs 
coupled to superconducting resonator. 



V. IMPERFECTIONS IN IMPLEMENTATION 

Here, we discuss different types of imperfections which 
can arise during the implementation of Grover's Algo- 
rithm. The relevant parameters during the implemen- 
tation are coupling constant gi, decay rate Fs of level 
|3) and cavity decay rate k. We consider spontaneous 
decay from the intermediate level |3), during resonant 
interaction of cavity field with |2) o |3) transition of 
SQUIDs (1) and (2). Here, we assume that under the 
condition that no photon from spontaneous emission is 
detected, conditional Hamiltonians for evaluation of sys- 
tem is given by [25|] 



1000) |0), 
1002) |0), 
1010) |1), 
1012) |1), 



1000) |0), 

1001) |0), 
1010) |1), 

-|oii)|i), 



1020) |1), 

1021) |1), 

1010) |1), 

1011) |1), 



%.at|2)^(3| + i/.c,)-iF3|3),(3|. (16) 



(13) 



Step 7. Apply microwave pulse (with 1^20^2 = 7''/2 
and phase ip — tt/2) to the SQUID 2 while a microwave 
pulse (with ^20^2 = 7r/2 and phase (p = — 7r/2) to the 
SQUID 1. The transformation |2)2 (10)2) -> |0)2 (- 12)2) 
for SQUID 2 while 12)^ (|0)i) ^ - (|2) J for SQUID 
1 is obtained. 

Step 8. Perform reverse operation of step 1. First wait 
for time interval ti — Tr/2gi during which resonator in- 
teracts resonantly to the |2) |3).^ transition of SQUID 
1. Then apply microwave pulse (with flsiti — 7^/2 and 
phase ip — tt) to SQUID 1. The overall step can easily 
be written as |2)j^ |1)^ |1)-^ |0)^. However, the state 
|0)j |0)^ remains unchanged. After applying steps 7-8, 
the system evolves 



The Hamiltonian given by Eq. (|16p is valid as long as 
the decay rate of level |2) is much smaller than F3. The 
decay rate T3 of SQUID (1) and (2) can effect the perfor- 
mance of three-qubit phase gate. Suppose each SQUID 
is initially prepare in generic state cos:^|0) -I- sini/|l) . 
Now the state of three SQUID qubit becomes \tp{0)) = 
a |000) + b |001) -I- c |010) + d |011) + e |100) + / |101) + 
g 1 110) + h\lll) , where the coefficients a — cos'^ ly, b = 
c = e = cos^ v sin v, d = f = g = cos u sit? v and 
h — sin^ 1^, satisfy the normalization condition. If we con- 
sider F3 = 0, then the state of the system after phase gate 
operation becomes |V'id(T)) = a |000) + b |001) + c |010) + 
rf|011)-f e|100) + /|101)+.g|110)-/i|lll). However, ff 
the decay of level |3) is included during phase gate then 
the expression for \4^{t)) becomes rather complex, there- 
fore it is not reproduced here. Fidelity of three-qubit 
phase gate is given by 
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FIG. 3: Average fidelity Fave of three-qubit controlled phase 
gate as a function of coupling constant Fg/ij. 



F 



(17) 



= 1 1 + cos^ V sin^ iy(r^ sin^ vir"^ -1)1^(15 



Here, r = {2gi/X)e~^sin9 with A = ^^igf - r|, rj = 
TrTs/Agi, and 9 = TrX/Agi. In order to realize the effects 
of decay on the performance of three-qubit phase gate, 
average fidelity over all possible three-qubit initial states 
is calculated which is given by 
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After carrying out integration, we obtain 



FIG. 4: Simulations results for three-qubit Grover's algorithm 

(19) 

involving spontaneous decay rate Fs from intermediate level 
|3). (a) Probability of finding the target state in case of 

Ta/g = 0, 0.001, and 0.004. (b) Fidelity of the state searched 
for F3/5 = 0.001 and 0.004. 



Fave = — (63 -I- 48r^ + 164r'* + 32r'^ + 8r^). (20) 
315 

Next, we show the plots of average fidelity as a function 
ofr^/g in Fig. [31 It can easily be verified that for = 
0, one has r = 1, which leads to F = Fave = 1. It is 
clear from Fig. [3] that average fidelity decreases due to 
the increase in the cavity decay rate. We also calculate 
the success probability of three-qubit phase gate which 
is given by 

P = (|a|' + |6|'-H|c|V|d|Vr*|e|' 

+r'\f\' + r'\9f + r'\h\'). (21) 

If we consider i' — 7r/4, then the corresponding success 
probability of three-qubit phase gate reduces to P = (4 + 

We take into account success rate of three-qubit phase 
gate (Eq. BT|) and Grover search itself. We consider 
<7i 52 ~ 5j without loss of generality and perform sim- 
ulation for finding the target state for T^/g = (ideal 



case), 0.001, and 0.004. Here, we consider the success 
probability for the case of target state |111) as shown in 
Fig. mja). It is clear from Fig. |4l[a) that success rate 
becomes closer to the ideal case when decay rate is suffi- 
ciently smaller than coupling constant i.e., T^/g = 0.001. 
Typical decoherence rate of SQUID is of the order 10^s~^ 
to 0.2 X 10^s~^ [ij, [l^ while coup ling constant can be 
achieved upto the value 3 x 10^s~^ [23|- This shows that 
for these parameters, three-qubit controlled phase gate 
and Grover's iterations can be performed with high fi- 
delity. Probability of success is highest at 6th iteration 
for ideal case. However, we should prefer 2nd iteration 
in the presence of dissipation because it has the highest 
value of fidelity. The effect of the level decay on the fi- 
delity of search state during iteration is shown in Fig. 
lU^b). It is clear from the Fig. IH^b) that fidelity of state 
to be searched decreases much rapidly for the case of 
larger decay rate as compared to smaller decay rate as 
the number of iterations increase. 

Next we consider the effect of cavity decay during the 
implementation of Grover's algorithm. If we assume that 
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the Q factor of superconducting cavity is infinite. Then, 
during the implementation of three-qubit phase gate in 
Sec. IV, transition of SQUID (1) from level |3) to |2) 
would result in the emission of one photon (See step 1 in 
Sec. IV). Transition of SQUID (2) from level |2) to |3) 
would absorb this photon with unit probability in step 3 
and vice versa for step 6 and 8. As a result, occupation 
probability of level |2) and |3) of SQUID (1) and (2) 
should be exactly one. However, if cavity relaxation is 
taken into account, then the occupation probabilities are 
expected to decay exponentially. Under the assumption 
that no photon actually leaks out during implementation 
time, we can write the conditional Hamiltonian as |8| 



h{gM^ |2) . (3| + H.c, ) - iKfll'a, 



(22) 



where, k is the cavity decay rate and gi is the coupling 
constant. For gi ^ k, we only need to focus on time evo- 
lution of the system governed by conditional Hamiltonian 
Eq. (I22p under the assumption of strong coupling limit. 
The fidelities and corresponding success probabilities of 
three-qubit phase gate can easily be obtained which are 
given by 



(4 + 2e~^'^* + e~'^* + e~3Kt/2-^2 
8(4 + 2e-4'^t + + e 



and 



P= 1(4 + 26-'^''* + 6-^''* + e 



(23) 
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FIG. 5: Simulation results for three-qubit Grover's algorithm 
involving cavity decay rate k (a) Probability of finding the 
target state |111) in case of k = 0, n/g = 0.004, n/g — 0.007. 
(b) Fidelity of the state searched for case n — 0, n/g — 0.004, 
K/g = 0.007. 



For K = 0, we have Fave = P = I - The success probability 
of target state |111) for ideal case and in the presence of 
cavity decay is shown in Fig. [Sja). It can be seen that 
the probability of finding target state decreases in the 
presence of cavity decay rate. The effect of cavity decay 
rate on fidelity of Grover's search iteration is also shown 
in Fig. Mjo)- It is clear from Fig. [S^b) that the fidelity 
of Grover's search decreases as a function of iterations 
much rapidly for higher cavity decay rate as compared 
to the lower cavity decay rate. 

It is clear from Fig. |4] and Fig. [5] that our proposed 
Grover's iterations can be performed with high success 
probability and fidelity as long as cavity decay rate and 
decay rate of the intermediate level |3) is small enough. 
The typical values of cavity decay rate is ~ 0.76/is 
{Q ~ 10^) . It may be mentioned that more rigorous 
analysis is required for the case of very low Q resonators. 



VI. DISCUSSION AND CONCLUSION 

We have discussed the implementation of 3-qubit 
Grover's algorithm using 4-level SQUIDs subjected to 
quantized and classical microwave fields. Grover's algo- 
rithm involves three-qubit phase gates and single-qubit 



gates. Here, we briefly estimate total operational time 
for three-qubit controlled phase gate. The total imple- 
mentation time is the sum of all interaction times in- 
volved in three-qubit controlled phase gate operation i.e.. 



T = 2ti + 2ti + 2t2 + 2t^ 



2U 



2t4 + t/. On substitut- 



ing the values of interaction times given in Sec. IV we 
obtain. 



1 



1 



1 



1 1 

+ ■ 



1 



<- ^ (25) 

^13 91 ^02 92 i^30 ^12 93 



Here, we consider without loss of generality ~ (72 ~ 
(73 ~ 3 X 10^ which is the same as given in Ref p7l |. 
Choosing = lOga, ^02 ~ O13 ~ f2i2 ~ ^30 ~ IO51, 
the operational time for three-qubit phase gate comes 
out to be T ~ 13ns. The operation time for single-qubit 
gate is about 1.5ns which can be applied to each qubit, 
simultaneously. Therefore, the estimated time for the im- 
plementation of three-qubit Grover algorithm performing 
two iterations is approximately 66ns. 

We studied the efffect of decay of level |3) using Eq. 
which is valid as long as the decay rate of level |2) 
(r2) is much smaller than T3. For example, T'^^ ~ 40/j,s 
and ~ 100/is, i.e., relaxation time of level |2) is larger 
than level |3) [26j. Moreover, during the phase gate op- 
erations given in Sec. IV, level |2) of SQUID 1 and 2 is 
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occupied through apphcation of SQUID-pulse resonant 
interaction and SQUID-resonator resonant interaction. 
Operation time t2 + t[ for SQUID 1 and i2+t3 for SQUID 
2, in either step is equal to 7r/(2r22o) + T^/{'^9i), which 
can be shortened by increasing the Rabi frequency r220 
and coupling constant gi. For the choice of ^20 = lOgi, 
and 5i ~ 3 X 10^s~^ as given in Ref [l^l, we have 
^2 + ^1 ~ ^2+^3 ~ 0.6ns, which is much shorter than 
Secondly SQUID can be designed to have long re- 
laxation time for level |2) , such that Fj"^ ^ t2 + ti. Thus 
decoherence due to relaxation of level |2) of SQUID 1 and 
2 is negligibly small. In our scheme direct coupling be- 
tween level |1) and |0) is not needed. Potential barrier 
between level |1) and |0) can be adiusted such that de- 
cay of level |1) is negligibly small [13 • Therefore, storage 
time of each qubit can also be made much longer. 

When levels |2) and |3) are manipulated by microwave 
pulses, resonant interaction between cavity field mode 
and 1 2) o |3) transition of each control SQUID is un- 
wanted. This effect can be minimized by setting the con- 
dition » g for SQUIDs (1) and (2). 

During the application of microwave pulse, off-resonant 
interaction of cavity field with |2) -H- |3) transition, in the 
presence of single photon inside the cavity is unwanted. 
It induces an unwanted phase which can effect the perfor- 
mance of the desired gate. This effect can be negligible 
under the condition V,i2 » gl/Ac for SQUID 3. How- 
ever, if this effect is included during the steps 4 and 5 for 
three-qubit phase gate, then the corresponding fidelity 
can easily be obtained. Here, we do not present the ex- 
pression for \iP{t)) owing to its complexity. The fidelity 
is given by 

F{x) = \l-x-x{a^-/3^-j^))\^, (26) 

where, x = a = cos^, /3 = S sin^/ 

7 = 2rJi2sin^/^(52 with S = gf/Ac and £, = 

TT-y/J^ + 4ri^2/4f^i2- Average fidelity over all possible ini- 
tial states is given by 

1 

Fave = J F{x)dx = ^{l + a^ + f3^+j^+j^ 


+f3^{l + 2j^)-a'^{l + 2(3'^ + 2j^)). (27) 

If off-resonant interaction during step 4 and 5 is not con- 
sidered, then we have (5 = 0,a = /3 = 0, and 7 = 1 which 
leads to F = Fave = 1- The plot of the average fidelity 
as a function of Rabi frequency ^12/93 is shown in Fig. 
[SI Wc choose Ac = 10(?3 for this plot. It can be seen 
that the average fidelity increases as a function of Rabi 
frequency ili2 applied to target SQUID. It is clear from 
Fig. [6|that for fli2 = 0.6(?3, we have Fave ~ 1- 

In our scheme, three qubit phase gate is implemented 
through interaction of the cavity field sequentially with 
|2) o |3) transition of selected SQUID either resonantly 



1 .000 r ' ^ ' n 

0.998 ; / ; 

. 0.996 / ^ 

0.994 / : 

0.992 ; ; 

0.990 ' ' ' 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 6: Average fidelity Fave of three-qubit controlled phase 
gate as function of Rabi frequency f2i2 for the case of Ac = 
Wgt. 

or off-resonantly as shown in step 1, 3, 4, 6, and 8 in 
section IV. It might be possible the transition channel 
between logical state |0) and auxiliary state |2) of differ- 
ent rf-SQUIDs could still be coupled by simultaneously 
interacting with common cavity field dispersively. These 
additional interactions may also induce the imperfections 
on desirable quantum operations. To study such kind of 
imperfections, we consider the case in which when cavity 
field interacts dispersively with |2) ^ |3) transition of 
SQUID 3 during the step 4 (wanted dispersive interac- 
tion), at that moment SQUID 1 and 2 arc coupled disper- 
sively with common cavity field via |0) O |2) transition 
and SQUID 3 with common cavity field via |1) o |2) 
transition (unwanted dispersive interactions) , simulta- 
neously. In this case, average fidelity of three qubit phase 
gate is given by 

Fave = —!— (50291 -I- 25960 COS Z + 
135135 

42778 cos(2Z) -I- 13888 cos(3Z) 
-F2208cos(4Z)), (28) 

here, Z = 7r(5Q2/Ao2)/(5|/Ac). Where, A02 is the de- 
tuning between |0) ^ |2) transition of SQUID 1 and 
2, and go2 is the coupling constant between the cavity 
field and |0) O |2) transition of the SQUID 1 and 2. For 
SQUID 3 we assume that A12 = A02, and f/12 = 302- 
If the unwanted off-resonant interactions are not consid- 
ered, then we have ggj/A^ = 0, Z = 0, which leads 
to Fave = 1- Fig- [7| shows the plot of average fidelity 
as function of Z. It is clear from Fig. [7| that fidelity 
decreases with increase in value of Z i.e., when the ef- 
fect of unwanted dispersive interactions starts to domi- 
nant on wanted dispersive interactions. To obtained a 
good fidelity, condition (702 /A02 ^ ffl/Ac should be ful- 
filled. For the choice of .go2/Ao2 — 0.06g|/Ac, we have 
Fave ~ 0.95. In principle, these requirements can be 
achieved through adjustment of level spacing of SQUIDs 
by varying the external fiux (jj^ or the critical current 
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FIG. 7: Average fidelity Fave of tliree-qubit controlled 
phase gate as function of Z — 7r((7o2/^02)/(g3/Ac). 



Ic (e.g., for variable barrier rf-SQUID) 17]. However, it 
may be mentioned that recently L. F. Wei et all., [2^ 
presented an alternative approach to robustly selective 
the interactions between the qubits and the cavity. 

Occupation probability of level |3) of SQUID 3 in step 4 
is P3 ~ Agl I (4(jr| + A^) [201 , needed to be negligibly small 
in order to reduce the gate error. We have ~ 0.04 for 
the choice of = lO^s which can be further reduce by 
increasing the ratio of Ac/ga. 

Here, we would like to mention that the physical im- 
plementation of three-qubit Grover's algorithm in cavity 
quantum electrodynamics (QED) has been proposed in 
a recent study by Yang et al @. The scheme is based 
upon the resonant interaction of three Rydberg atoms ini- 
tially prepared in a coherent superposition state travers- 
ing through a single-mode microwave cavity. As com- 
pared to the flying qubits, here we have considered sta- 
tionary qubits defined through two lowest level of four 
level SQUIDs. Our scheme is based on the generation 



and absorption of single photon in high Q cavity using 
SQUIDs. The generation of single microwave photon in 
superconducting qubit have been reported, experimen- 
tally in some recent studies et al. [30|. 

When slowly changing Rabi frequencies are applied to 
satisfy the adiabatic passage, gate times becomes slow 
i.e., gate time is on the order of 1 ms to a few microsec- 
onds [sij . However, in our scheme we do not required 
slowly changing the Rabi frequency, thus gate is signifi- 
cantly faster i.e., gate time is on the order of 13ns. How- 
ever, it may be mentioned that fast pulses may intro- 
duce new imperfections, for example, precisely designing 
of the duration of pulses might not be easy. However, 
adiabatic is not always slow. Recently, an interesting 
proposal based on Stark-chirped rapid adiabatic passage 
to resolve the issue of precise pulse duration have been 
proposed [32j . 

In conclusion, we have proposed a scheme for the real- 
ization of three-qubit Grover's algorithm based on three- 
qubit phase gate using four-level SQUIDs coupled to 
single-mode superconducting resonator. In this proposal, 
adjustment of level spacing during the operation, slowly 
changing Rabi frequencies (to satisfy adiabatic passage) , 
and the use of second-order detuning (to achieve off- 
resonance Raman coupling between two relevant levels) 
are not required. Thus, implementation time is signif- 
icantly faster and has operation time of the order of 
nanoseconds. The coupling constants of each SQUID 
with the resonator are different. Thus, an unavoidable 
non-uniformity in device parameters can be accommo- 
dated. We consider the effect of imperfections in the 
system which include the decay of the cavity field and 
the relevant level. We also incorporate the influence of 
unwanted off-resonant interaction during the gate imple- 
mentation. Our results show that the marked state can 
be searched with high fidelity even in the presence of im- 
perfections in the system. 
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